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Abstract 

We present measurements of time-dependent CP- violation asymmetries for the decays B° — > tj'K . 
The data sample corresponds to 347 million BB pairs produced by e + e~ annihilation at the T(4S) 
resonance in the PEP-II collider, and collected with the BABAR detector. The preliminary results 
are S = 0.55 ± 0.11 ± 0.02, and C = —0.15 ± 0.07 ± 0.03, where the first error quoted is statistical, 
the second systematic. 
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1 INTRODUCTION 



Measurements of time-dependent CP asymmetries in B° meson decays through a dominant Cabibbo- 
Kobayashi-Maskawa (CKM) favored b — > ccs amplitude [1] have provided a crucial test of the mech- 
anism of CP violation in the Standard Model (SM) [2]. For such decays the interference between 
this amplitude and B°B° mixing is dominated by the single phase (3 = arg (—V c dV* b /VtdV t l) of the 
CKM mixing matrix. Decays of B° mesons to charmless hadronic final states such as 7]'K° proceed 
mostly via a single loop (penguin) amplitude with the same weak phase as the b — > ccs transition 
[3], but CKM-suppressed amplitudes and multiple particles in the loop introduce additional weak 
phases whose contribution may not be negligible [4, 5, 6, 7, 8]. 

For the decay B° — > rj'K , these additional contributions are expected to be small within 
the SM, so the time-dependent asymmetry measurement for this decay provides an approximate 
measurement of sin2/3. Theoretical bounds for the small deviation AS between the time-dependent 
CP-violation parameter S measured in this decay and in the charmonium-i^° decays have been 
calculated with an SU(3) analysis [4, 5] from measurements of B° decays to pairs of neutral light 
pseudoscalar mesons [9, 10]. The most stringent of these is given by Eq. 19 in [5], which assumes 
negligible contributions from exchange and penguin annihilation, and has a theoretical uncertainty 
less than ~0.03. With newer measurements [11] we obtain an improved bound AS < 0.08 [12]. 
QCD factorization calculations conclude that AS is even smaller [7]. A significantly larger AS 
could arise from non-SM amplitudes [8]. 

The time-dependent CP-violation asymmetry in the decay B° — ► rj'K has been measured 
previously by the BABAR [13] and Belle [14, 15] experiments. In this paper we update our previous 
measurements with an improved analysis and a data sample 1.5 times larger. 

2 THE BABAR DETECTOR AND DATASET 

The data were collected with the BABAR detector [16] at the PEP-II asymmetric-energy e + e~ 
collider [17]. An integrated luminosity of 316 fb -1 , corresponding to 347 million BB pairs, was 
recorded at the T(4S) resonance (center-of-mass energy y/s = 10.58 GeV). 

Charged particles from e + e~ interactions are detected, and their momenta measured, by a 
combination of five layers of double-sided silicon microstrip detectors and a 40-layer drift chamber, 
both operating in the 1.5 T magnetic field of a superconducting solenoid. Photons and electrons are 
identified with a CsI(Tl) electromagnetic calorimeter (EMC). Further charged particle identification 
(PID) is provided by the average energy loss (dE/dx) in the tracking devices and by an internally 
reflecting ring imaging Cherenkov detector (DIRC) covering the central region. The instrumented 
flux return (IFR) of the magnet allows discrimination of muons from pions. 

3 ANALYSIS METHOD 

3.1 Time evolution of a B°B° pair 

From a candidate BB pair we reconstruct a B° decaying into the CP eigenstate / = rj'K® or 
/ = rj'K® (Bcp). From the remaining particles in the event we also reconstruct the vertex of the 
other B meson (Stag) an d identify its flavor. The difference At = top — t tag of the proper decay 
times tcp and ttag of the signal and tag B mesons, respectively, is obtained from the measured 
distance between the Bqp and B tag decay vertices and from the boost (/?7 = 0.56) of the e + e~ 
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Table 1: Selection requirements on the invariant masses of resonances and the laboratory energies 
of photons from their decay. 



State Invariant mass (MeV) £'(7) (MeV) 

120 < m(77) < 150 > 30 

120 < 771(77) < 155 > 30 

490 < 771(77) < 600 > 50 

520 < 777(7r+7T-vr ) < 570 
945 < m(Tr + TT-f]) < 970 
930 < T7i(7r+7r-7) < 980 > 100 
470 < 77i(vr+7r-) < 980 
486 < T7i(7r+7r-) < 510 
468 < 77i(7r°7r°) < 528 



system. The At distribution is given by: 

e -|At|/r 

F(At) = — [1 t Au; ± (1 - 2w) (Ssm(Am d At) - Ccos(Am d At))]. (1) 

4r 

The upper (lower) sign denotes a decay accompanied by a B° (B°) tag, r is the mean B° lifetime, 
Am^ is the mixing frequency, and the mistag parameters w and Aw are the average and difference, 
respectively, of the probabilities that a true B° is incorrectly tagged as a B° or vice versa. The 
tagging algorithm [18] has six mutually exclusive tagging categories based on quantities such as the 
sign of charge of a lepton, kaon, or soft pion from D*, grouped according to their response purities. 
The measured analyzing power, defined as efficiency times (1 — 2w) 2 summed over all categories, is 
(30.4±0.3)%, as determined from a large sample of B-decays to fully reconstructed flavor eigenstates 
(-Sflav)- The parameter C measures direct CP violation. If C = 0, then S = — 7)sin2/3 + AS, where 
r] is the CP eigenvalue of the final state (—1 for rj'K®, +1 for rj'K®). 

3.2 Event selection 

We establish the event selection criteria with the aid of a detailed Monte Carlo (MC) simulation 
of the B production and decay sequences, and of the detector response [19]. These criteria are 
designed to retain signal events with high efficiency. Applied to the data, they result in a sample 
much larger than the expected signal, but with well characterized backgrounds. We extract the 
signal yields from this sample with a maximum likelihood (ML) fit. 

The i?-daughter candidates are reconstructed through their decays ir° —>■ 77, 77 — > 77 (t/ 77 ), 
7/ -► vr+vr-vr (rj 3n ), rf -► 7/ 77 7r+7r~ (r^ (ry)7r7r ), rf -► m^ir" (^ (3w)ww ), rf -► p°7 (r/„), where 

p° — > 7T + 7r _ , Kg — > 7r + 7r _ (i^gj ) or 7r°7r° (K® 00 ). Table 1 lists the requirements on the invariant 

mass of these particles' final states. Secondary charged pions in 7/ and rj candidates are rejected 
if classified as protons, kaons, or electrons by their DIRC, dE/dx, and EMC PID signatures. We 
require candidates to have a flight length with significance >3a. Signal K® candidates are 
reconstructed from clusters of energy deposited in the EMC or from hits in the IFR not associated 
with any charged track in the event. From the cluster centroid and the B° decay vertex we 
determine the direction (but not the magnitude) of the K\ momentum p x o . 



7r (from 77 37r ) 
vr (from K° s00 ) 



VriTTTT 

K° 
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For decays with a K® we reconstruct the S-meson candidate by combining the four-momenta 
of the Kg and 7/ and imposing a vertex constraint. Since the natural widths of the 77, rf, and 
7T° are much smaller than the resolution, we also constrain their masses to world-average values 
[20] in the fit of the B candidate. From the kinematics of T(4S) decay we determine the energy- 
substituted mass toes = \J {\s + po • pb) 2 /Eq — p B and the energy difference AE = E* B — \\fs, 
where (£b,Po) and (Eb,Pb) are four-momenta of the T(4S) and the B candidate, respectively, 
and the asterisk denotes the T(4S) rest frame. The resolution in toes is 3.0 MeV and in AE is 
20 — 50 MeV, depending on the decay mode. We require 5.25 < toes < 5.29 GeV and \AE\ < 0.2 
GeV (-0.01 < AE < 0.08 GeV for B° -► r/K®). 

For a B° — > rfK® candidate we obtain AE and p K o from a fit with B° and K® masses con- 
strained to their accepted values [20]. To make a match with the measured K® direction we 
construct the missing momentum p m i ss from po and all charged tracks and neutral clusters other 
than the K\. We then project p m iss onto P^o, and require the component perpendicular to the 

beam line, P^issJ.' 

to satisfy p^ ± -p K o ± > -0.5 GeV/c. This value reflects the resolution, and 
was chosen to minimize the yield uncertainty in the presence of continuum background. 

For all B° — > r/'K candidates we require for At and its error a^t, \At\ < 20 ps and a^t < 2.5 

ps. 

3.3 Background rejection 

Backgrounds arise primarily from random combinations of particles in continuum e + e~ — ► qq events 
(q = u, d, s, c). We reduce these with requirements on the angle 9t between the thrust axis of the 
B candidate in the T(4S) frame and that of the rest of the charged tracks and neutral calorimeter 
clusters in the event. The distribution is sharply peaked near | cos 6t\ = 1 for qq jet pairs, and nearly 
uniform for i?-meson decays. The requirement, which optimizes the expected signal yield relative 
to its background-dominated statistical error, is | cos#t| < 0.9 (| cos#t| < 0.8 for B° — > rj'K®). 

In the ML fit we discriminate further against qq background with a Fisher discriminant T 
that combines several variables which characterize the production dynamics and energy flow in the 
event [10]. It provides about one standard deviation of separation between B decay events and 
combinatorial background. 

For the rj^ decays we require |cos^ cc | < 0.9 to exclude the most asymmetric decays where 
soft-particle backgrounds concentrate and the acceptance changes rapidly. Here Q p Acc is the angle 
between the momenta of the p° daughter 7r~ and the rj', measured in the p° rest frame. 

For B° — ► rfK® candidates we require that the cosine of the polar angle of the total missing 
momentum in the laboratory system be less than 0.95, to reject very forward qq jets. The purity of 
the K\ candidates reconstructed in the EMC is further improved by a requirement on the output 
of a neural network (NN) that takes cluster-shape variables as its inputs. The NN was trained 
on MC signal events and data events in the sideband 0.04 < AE < 0.08 GeV. We checked the 
performance of the NN with candidates in the larger B — > J/tpK® sample. 

The average number of candidates found per selected event is in the range 1.08 to 1.32, de- 
pending on the final state. We choose the candidate with the smallest value of a x 2 constructed 
from the deviations from expected values of one or more of the daughter resonance masses, or with 
the best vertex probability for the B, depending on the decay channel. In B° — > rj'K® if several 
B candidates have the same vertex probability, we chose the candidate with the K\ reconstructed 
from, in order, EMC and IFR, EMC only, or IFR only. From the simulation we find that this 
algorithm selects the correct-combination candidate in about two thirds of the events containing 
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multiple candidates, and that it induces negligible bias. 
3.4 Maximum likelihood fit 

We obtain the common CP-violation parameters and yields for each channel from a maximum 
likelihood fit with the input observables AE, tbeSj an d At. The selected sample sizes are given 
in the first column of Table 2. Besides the signal events they contain qq (dominant) and bb with 
b — > c combinatorial background, and a fraction that we estimate from the simulation to be less than 
1.1% of cross feed from other charmless BB modes. The charmless events (henceforth refered to as 
BB) have ultimate final states different from the signal, but similar kinematics, and exhibit broad 
peaks in the signal regions of some observables. We account for these with a separate component in 
the probability density function (PDF). For each component j (signal, qq combinatorial background, 
or BB background) and tagging category c, we define a total probability density function for event 
i as 

Pj >c ee • VjiAE*) • V^r) ■ V 3 (At\ a\ t - c) , (2) 

except for B° — > rj'K® for which Vj(m^ 1 $ i ) is omitted. The factored form of the PDF is a good ap- 
proximation, particularly for the combinatorial qq component, since correlations among observables 
measured in the data (in which qq dominates) are small. Distortions of the fit results caused by 
this approximation are measured in simulation and included in the bias corrections and systematic 
errors discussed below. 

With Yj defined to be the yield of events of component j, and / J)C the fraction of events of 
component j for each category c, we write the extended likelihood function for all events belonging 
to category c as 

N c 

£ c = exp ( - ]T Yjhc) U(W-,<P: i; ,, + Y m f q ^Vi q + Y B b!bb, c V B b\ (3) 

3 i 

where N c is the number of events of category c in the sample. We found that the BB background 
component is needed only for the channels with rj^. We fix both / s i g , c and /bb,c t° /B flav ,c 5 the 
values measured with the large -Bfj av sample [21]. The total likelihood function Cd for decay mode 
d is given as the product over the six tagging categories. Finally, when combining decay modes we 
form the grand likelihood C = Y\ C-d- 

The PDF "P s i g (At, 0"At,c) is given by F(At) (Eq. 1) with tag category (c) dependent mistag 
parameters convolved with the signal resolution function (a sum of three Gaussians) determined 
from the -Bfl av sample. We determine the remaining PDFs for the signal and BB background 
components from fits to MC data, for which the resolutions in AE and tties ar e calibrated with 
large control samples of B decays to charmed final states of similar topology (e.g. B — > D{Kttk)tt). 
For the combinatorial background the PDFs are determined in the fits to the data. However we 
first deduce the functional form from a fit of each component alone to a sideband in (mES; AE 1 ), so 
that we can validate the fit before applying it to data containing the signal. 

These PDF forms are: the sum of two Gaussians for 7 3 s ig(^Es) and V S \ Z {AE); the sum of three 
Gaussians for V q q(At; c); a conjunction of two Gaussian segments below and above the peak with 
different widths for Vj(T) (a small "tail" Gaussian is added for V q q(!F))\ a linear dependence for 
Vqq(AE); and for V q q(mEs) the function xV 1 — x 2 exp [— £(1 — x 2 )], with x = 2mEs/v / ^- These 
are discussed in more detail in [10]. 

We allow the parameters most important for the determination of the combinatorial background 
PDFs to vary in the fit. Thus for the six channels listed in Table 2 we perform a single fit with 
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109 free parameters: —rjS, C, signal yields (6), rj' p ^K° BB background yields (2), continuum 
background yields (6) and fractions (30), background At, ities, AE, J- PDF parameters (63). The 
parameters r and Am,; are fixed to world-average values [20]. The symbol S refers to S v , K o, and 
inclusion of the CP eigenvalue ry of the final state accounts for the expected difference in sign with 
respect to rfK®. 

We test and calibrate the fitting procedure by applying it to ensembles of simulated qq experi- 
ments drawn from the PDF into which we have embedded the expected number of signal and BB 
background events randomly extracted from the fully simulated MC samples. We find negligible 
bias for C. For S we find and apply multiplicative correction factors for bias from dilution due to 
BB background, equal to 1.02 in the final states rj' frj K^ + ^_ and rf^ W7r K®, and 1.05 in J^-K^o^o- 

3.5 Fit Results 



Table 2: Results with statistical errors for the B° — ► rj'K time-dependent fits. 

Mode Events to fit Signal yield — rjS C 

■tLZa K" + _ 612 206 ± 16 0.60 ± 0.24 -0.26 ± 0.14 

ri'K , _ 10905 503 ±28 0.50 ±0.15 -0.26 ±0.11 

'Pi 7T^7T 

V'^^K^ 164 63 ±8 0.85 ±0.38 0.24 ±0.26 

<(ry)™*2%° 446 50 ± 9 0.77 ±0.44 -0.25 ± 0.36 

n'K% o 12559 114 ± 23 0.42 ± 0.47 0.30 ± 0.30 









0.57 ±0.11 


-0.18 ±0.08 




3389 


168 ± 21 


0.39 ±0.30 


0.20 ±0.23 


rj'K" 






0.55 ±0.11 


-0.15 ±0.07 



CM 

o 

> 

CM 



200- 



HlOOh 
a> 

5 o 



(a) BABAR 

preliminary 




5.25 5.26 5.27 5.28 5.29 
m ES (GeV/c 2 ) 



> 

0) 



o 

CM 
(0 



200- 



g 100h 
> 

LU 



xBaBar 

preliminary 




.2 -0.1 0.1 0.2 
AE (GeV) 



Figure 1: Distributions projected onto (a) thes and (b) AE for B — > r]'Kg candidates 



Results from the fit for the signal yields and the CP parameters S and C are presented in Table 
2. In Fig. 1 we show for B° — > r]'K® the projections onto tues an d AE for a subset of the data 
for which the signal likelihood (computed without the variable plotted) exceeds a mode-dependent 
threshold that optimizes the sensitivity; the corresponding distribution in AE for B° — > vj 
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Figure 2: Distribution projected onto AE for B° — > r/'K® candidates. Points with error bars 
represent the data, the solid line the fit function, and the dashed line its background component. 



given in Fig. 2. Fig. 3 gives the At projections and asymmetry of the combined modes for events 
selected as for Figs. 1 and 2. We measure a correlation of 3.0% between S and C in the fit. 

We perform numerous crosschecks of our fitter: time-dependent fits for B + decays to the charged 
final states rj' rj ^^ 7T7T K + , rj' p ^K + , and rj'^^^K^ ; fits removing one fit variable at a time; fits without 
BB PDFs; fits with multiple BB components; fits allowing for non-zero CP information in BB 
events; fits with C = and others. In all cases, we find results consistent with expectation. The 
value S^ftO = 0.57 ±0.11 is larger than our previous measurement S v , K o = 0.30 ± 0.14 [13] as 
a result of the larger data sample and events added or removed as a result of changes in the 
reconstruction and selection. For events common to the two datasets we find close agreement of 
the values of S and C. 
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error bars), fit function (solid line), and background function (dashed line), for (a, d) B° and (b, 
e) B° tagged events, and (c, f) the asymmetry between B° and B° tags. 
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3.6 Systematic studies 

We find systematic uncertainties from several sources (in decreasing order of magnitude) : variation 
of the signal PDF shape parameters within their errors, modeling of the signal At distribution, use 
of At signal parameters from the i?fl av sample, interference between the CKM-suppressed b — ► ucd 
amplitude and the favored b — ► cud amplitude for some tag-side B decays [22], BB background, 
SVT alignment, and position and size of the beam spot. The -Bfl av sample is used to determine 
the errors associated with the signal At resolutions, tagging efficiencies, and mistag rates. We take 
the uncertainties in tb and Anid from the published measurements [20]. Summing all systematic 
errors in quadrature, we obtain 0.02 for S and 0.03 for C. 

4 RESULTS AND DISCUSSION 

In conclusion, we have used samples of about 940 B° — > rj'K® and 170 B° — > r/'K® events to 
measure the time-dependent CP violation parameters in B° — > r/'K S = 0.55 ± 0.11 ± 0.02 and 
C = —0.15 ± 0.07 ± 0.03. Our result for S is consistent with the world average of those measured 
in B° — > J/t/jKg [18, 15], and inconsistent with zero (CP conservation) by 4.9 standard deviations. 
Our result for the direct-CP parameter C is 1.8 standard deviations from zero. The results are 
preliminary. 
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